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1 INTRODUCTION Theoretical studies of energetic particles in 
space are often referred to a special frame of reference. To compare 
theory with experiment, one has to transform the particle distribution 
from the special frame to the observer’s frame, or vice versa. Various 
methods have been derived to obtain the directional distribution in the 
comoving frame from the directional fluxes measured on a spacecraft 
(e.g., Erdos, 1981; Sanderson et al., 1985, and references therein). 
These methods have become progressively complicated as increasingly 
detailed directional particle data become available. 

We present here a set of 2nd order correct formulae for the 
transformation of an arbitrary differential intensity distribution, 
expressed as a series of spherical harmonics, between any two frames in 
constant relative motion. These formulae greatly simplify the 
complicated procedures currently in use for the determination of the 
differential intensity distribution in a comoving frame. 


2. THE DESIRED FORM OF THE TRANSFORMATION In the observer's frame, 
let the particle differential intensity w.r.t. momentum (or rigidity) be 
expressed in terms of the spherical harmonics CR, Sjj as : 
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We wish to transform this to the frame moving with velocity W relative to 
the observer , also m terms of spherical harmonics : 
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PP?C > = associated Legendre function, p = particle momentum, and 0, 0 are 
the polar and azimuthal angles of the particle arrival direction. The 
superscript or subscript S indicates that the quantity is referred to the 
frame of the observer /spacecraft. 

In short, our purpose is to find the relation between the coeffi- 
cients A„ ra , B tirn in the comoving frame and the coefficients A^ fe , B5„ 
and their p-denvati ves in the observer's frame, and vice versa. 

We denote below the components of W along and normal to the magnetic 
field by W| t and W ± respectively, and define Wy= W„ /W, W ± = W x /W^ and 
€ - W/v, with v the particle speed. We choose the polar coordinate axis 
to point along the magnetic field, and measure 0 from the direction of 
Wj_. Hence, »), 0 denote the particle pitch-angle and gyrophase respect- 
ively. We will refer to this coordinate system as the standard one. 


3* _ SOME CURRENT SOLUTIONS Gleeson and Ax ford (1968) and Forman 
(1970) obtained 1st order transformations from (2) to d), for j„=A 0 o(p). 
Forman’s elegant approach uses the Lorentz invariance of the phase-space 
distribution function. Balogh et al. (1973) obtained a 2nd order result 
for the isotropic case and a 1st order result for the anisotropic case 
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that includes the n=l terms. Ng (1984) obtained a 2nd order 
transformat ion for gyrotropic distr lbutions (j p = X A„o(p)FV%(cosQ) ) valid 
for relativistic particles. 1 

Another approach, taken by Sanderson et aim (1985), is based on a 
method by Ipavich (1974). This is outlined below. One may use eqn 
(12) of Ipavich (1974) to transform each term in the series (1), provided 
one pays the price of using a non-standard coord mate system in which the 
polar axis points along the velocity W. Denote the non- 
standard polar and azimuthal angles by 9*, 0*. This allows any term 
ASm<‘p©'*CE(9jg, 0g) for the phase-space density to be transformed to: 

A S lp(l-2ecos0 -t-e ) 5 ]C m {tan [sin9 /(cos0 -e),<)> ]}, (4) 
nm n 

and similarly for B® rft (p a )CfJ(9 §,0g). Following Ipavich and Sanderson 
et al. f one assumes a power-law spectrum (p Q d/dp e ) log A® n ,(p s ) = - 2r nf 
and obtains from (4), _ r -i 

(l-2ecos0 +e 2 ) n ^m^^ C n^ tnn ^ f sin0 /( cos0*-e ),((>*]} . ( 

Eqn t’5) may be compared to Sanderson et aim 9 s eqn (All). There remains 
the task of converting (5) to the form (2), including a change from 
9^,0^ back to the standard angular coordinates 9,0. 

Sanderson et al, accomplish the above by numerically integrating the 
transformed intensity obtained from (5) (their eqns (All) and (A12)) over 
each sector, each energy interval and each telescope to produce an 
equivalent set of counts m the comoving frame , and then fitting (2) to 
the equivalent set of counts. 

Although the assumption of a power-law spectrum may be avoided by 
using (4), the numerical integration and the second fitting are time- 
consuming and unnecessary for say, >10 KeV protons, because we are able 
to obtain the transformation posed in Section 2, correct to CKW/v)* and 
for arbitrary direction of W. 


4. DERIVATION OF THE TRANSFORMATION FORMULAE 


We begin with 


2 2 

j n (p.e,*> = ££ { *1 (p )c™(9 ,+ )+ b s ( P )s m (e ,♦ )}, 
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and consider each term in the series. First we transform the spherical 
harmonics OR (9 s ,0©). Galilean transformation of momentum gives 

cos 0^ = cos 0 - e(ft|| -t; cos 0) + (3C^-l)cos 0- 2W||C] + O(c^), (7) 


sin m 0 cos m<() = sin m 0 cos + me[ £ sin^O cos m(f) - W.. sin" x 0 cos (m-1 )<|> ] + 

S S _L 

2 rr , i . ttu . ~ 2 " _ . m- 1 . 


TO, 


.. TO-1 Q 


+ Je {[m(m+2)£ -m]sin 0 cos - 2m sin 0 cos(m-l)<f> + 


where 


+ W^(m-l)m sin m ^0 cos(m-2)<(>} + O(e^) , 

A A A A 

£ = - v.W = W|j cos 0+ W sin 0 cos <f> > 


t, cos m<|> = Wjj cos 0 cos m<|> + ^W^sin 0[cos(m-lH + cos(m+l)<J>] , 

2 /n 2 2 a 2 2 a a 

K cos m<|> = (Wjj cos 9 + sin 0) cos m<t> + sin 0 cos 0 x 
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x [cos(m l)<j> + cos(m+l)<fi] + ijW^ sin 2 0[cos(m-2 )<(> + cos(m+2)<|>] . (11) 

Using (7) - (101, and various properties of PR, we obtain 
C H ( W = C^(0,'t»)+e( n + 1 )U^ 1 ( e , < (, ) ( 0 , <(» ) + E 2 (n-l)(n+l)X™ 2 (e,<J>) + 


+ e n(n+l)Y n (0,<f>) + e n(n+2)Z^ 2 (0 ,<J>) + 0(e°), 

where (2n+l)U™ C = (n+m-1) (n+m)C m * - W (n+m)C m - Jw C™** 

n- J- l n-1 II n-1 2 1 n-1 
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(2n-l)(2n + l)X- = ^ - tf,, W, + 

+ * (2 "ll -”l ) T^SlT C n-2 + ^(n+n,)^ + £w 2 C™+ 2 , (15) 

(2n-l)(2n+3)Y™ C = £fl 2 {f-^2)"(n-m) ! °n 2 + 2m -DW,| ^(n-m+lHn+mJC™ -1 * 

- \ [(2n 2 +2n-2m 2 -l)W 2 + (n 2 +n+m 2 -l)W 2 ]c m +s(2m+l)W M W. C m+1 +^W 2 c ,n+ ^ IC'l 

11 1 n II 1 n M n* 
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+ «*K> C 2 - i»„ * &1CI • ,m 

Note that the functions U, V, X, Y f and Z are all expressed m terms of 
the spherical harmonics <O r 0) f and that 
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( 18) 
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C[(n-2)U'^ 1+ (n+3)^] - |c™ = (4-2n)X^ 2 -3Y™ C + (2n+6)Z™ 2 . < 20 ) 

Secondly, we expand A« m (p B ) about p = p a : 

A n» ( P« ) = A L(P)+^PA:!(P) + iE 2 [(l-C 2 )pA^(p) + C 2 p 2 A"!(p)]+0(e 3 ). (21) 


nm 


nm 


nm 


nm 


nm 


Then we multiply (21) and (12), and use (10), (11), (19) and (20) to 
express A£ m fp a )CR($ 0 , 0 S ) in terms of a series of CR<9,0) correct to 
0(£®)- This expression actually applies to the transformation of phase- 
space distribution function, but will not be exhibited here. 

Final ly, multiplying 


Vp /( VV = 1 + 3e ^ " (3e 2 /2)(l-5C 2 ) + 0(e 3 ) 
into the above expression, and simplifying as before, we obtain 

2 ,, 2 
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the right hand side being evaluated at (p,9,0). Eqns (8>, dO) - (17), 
(19) - (20) and (23) hold when we replace all the cosines dependent on 0 

by sines, CR by SR, UR2» by UR®! , etc. Substituting (23) into (6), we 
obtain the desired result in the form (2), as well as the relation 
between A„ m and Ag„ and their derivatives. These, however, are not 
exhibited here for lack of space. For the inverse transformation, we 
merely need to interchange quantities labelled and unlabelled with S, 
and reverse the signs of W|| and Wj_ in the above. 

S. SOME IMPLICATIONS For >10 KeV protons, once A% and their 1st 
and 2nd derivatives are experimentally determined, one need only substi- 
tute these into the formulae to obtain the relevant A„ ra in the comoving 
frame, bypassing the complicated procedures mentioned in Section 3. 

The formulae are useful in other regards. For example, we have 


A n0 = A n0 + e ^ n+1 )/^ 2n+3 ^Hn-l+pd/dp][iW 1 (n+2)A n+1 ^ 1 -W || ( n+1 > A n+1)0 l + 

+ e[n/(2n-l)][n+2-pd/dp][|W 1 (n-l)A nl ^ 1 +W|| nA nl ^ 0 ] + 0(e ). (24) 

For highly anisotropic solar particles, the gyrotropic terms A k0 » the 
non-qyrotropic term A*i. The former may be eliminated in *24.) by 
choosing the ExB drift frame (Hq ® 0). If in addition, An— i r i and An-*-i,i 
are neglible, then A n0 ~ AJJo = 0<& 2 ). 

Next consider the 1st harmonic terms* 


A 10 = A 10 + e ”ll (3A 00 JpA 00 _2pA 20 /5) " E ”l pA 21 + 0(e > » 

“ A ll = “ A ll + e ^l (3A 00 _pA 00 +pA 20 /5_6pA 22 /5) + 3e ”|| pA 21 /5 + ° (e2) ’ 
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" B ll = _B 11 + 3eW ||P B 2l /5 " 6eW iP B 22 /5 + 0(e ) * 

fen if all the non-qyrotropic terms vanish, if pA*»o f 0, 
order 1st harmonic anisotropy vector is not aligned with W 1 


t 25) 
the 1st 


6. REMARKS The new formulae are useful in the study of transverse 
anisotropies. The derivation in Section 4 is being extended to cover 
relativistic particles (see also Ng, 1984). These and other 
considerations and a more complete account of the present work will be 
presented elsewhere. 
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